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Q 1 ABSTRACT 
IT) ' 

■ We study the properties of the diffuse X-ray background by using the results of a cosmo- 

logical hydrodynamical simulation of the concordance ACDM model. The simulation follows 
gravitational and gas dynamics and includes a treatment of physical processes like radiative 
cooling, star formation and supernova feedback. From the simulation outputs, we produce 
a set of two-dimensional maps of the intergalactic medium X-ray emission integrated over 
redshift. We find that the signal in the soft (0.5-2 keV) band is lognormally distributed with 
£2 ' a mean intensity of about 4 x 10~ 12 erg s _1 cm~ 2 deg~ 2 ; approximately 40 per cent of the 

emission originates from warm-hot gas (defined as baryons with 10 5 < T < 10 7 K), and 90 per 
cent comes from structures at z < 0.9. Since the spectrum is soft, being mostly provided by the 
intergalactic medium at low temperature, the total mean intensity in the hard (2-10 keV) X-ray 
band is smaller by a factor of about 4. In order to constrain the physical processes included 
in our simulation, we compare our results with the observed upper limit (1.2 ±0.3) x 10~ 12 
erg s cm~ 2 deg~ 2 of the soft X-ray emission due to diffuse gas. To this purpose, we remove 
the contributions of observable extended objects (groups and clusters of galaxies) from the 
simulated maps by adopting different detectability criteria which are calibrated on the prop- 
erties of systems at intermediate redshifts observed by Chandra. We show that the simulated 
diffuse soft X-ray emission is consistent with the present observed upper limit. However, if 
future measurements will decrease the level of the unresolved X-ray background by a factor 
of two, a more efficient feedback mechanism should be required to suppress the soft emission 
of the gas residing within filaments and group-size haloes. 

Key words: diffuse radiation - intergalactic medium - hydrodynamics - galaxies: clusters: 
general - X-rays: general - cosmology - theory 



1 INTRODUCTION 

A general result of hydrodynamic simulations of cosmic structure 
formation is that, during the formation of galaxy clusters via grav- 
itational collapse, a large amount of gas remains out of these struc- 
tures at the present time. An important fraction of the total baryons 
of the universe (roughly 40-50 per cent) undergoes a process of 
shock-heating that starts at z ~ 2 and heats the baryons to an in- 
termediate temperature of 10 to 10 7 K: this gas phase is often re- 



ferred to as the Warm-Hot Intergalactic Medium (WHIM) (see, e.g., 
iDave et alJl200lt Icroft et alJl200lh . Simulations also show that 
these baryons are not uniformly distributed: they constitute a fila- 
mentary network linking the largest virialized objects, the so-called 
cosmic web ( Bo nd et alll99rJ) . The existence of these structures is 
believed to be the solution to the missing baryon problem (see, e.g., 
ICen & Ost riker 1999), namely the fact that more than half the nor- 
mal matter is yet to be detected by instruments: in fact, the density 
of filaments is expected to be very low and consequently their de- 
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tection is made difficult by both their extremely low surface bright- 
ness and projection effects (see, however, some first claimed detec- 
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tions in X-ray analyses: Zappacosta et al. 2002; Markevitch et al. 

l2003l l Finoguenov et al. 2003]). In any case, the intermediate tem- 
perature of the cosmic web gas suggests that it can produce a non- 
negligible X-ray emission through thermal bremsstrahlung. Given 
the physical state of the gas, this emission is expected to appear 
mostly as a diffuse background at soft energies. 

It is now generally accepted that some of the observed clus- 
ter X-ray properties, like the scaling relations between mass, X- 
ray luminosity and temperature and their redshift evolution, the 
temperature profiles, or the entropy excess in the central regions 
of poor clusters and groups, cannot be reproduced by a simple 
model of gasdyna mics based on gra vitational heating alone (see 
the discussion in iRosati et alJl2002t IVoTilEooll and references 
therein). In an attempt to explain these discrepancies, new pro- 
cesses have been included in the models of cluster formation, 
largely increasing the complexity of the physics of the intraclus- 
ter medium (ICM): radiative cooling, energy feedback from su- 
pernovae (SNa;) and active galactic nuclei (AGN), thermal con- 
duction, turbulence, magnetic fi elds, etc. (see, e.g., |Bowej[l997|; 
Cavaliere et alll998l:lBalogh et all 19991: IPearce et all2000l:lBrvarJ 
200filTozzi & Normarj20qitlMuanwong et alj200lUVoit & Brvanl 
2001bl: iBabul et alJl2002t IVoit et alj|2002t iTornatore et alj|200a 
Jubelsas et all2004lDolag et all2005l:lDi Matteo et all2005t) . All 
these processes affect the properties of both the gas within viri- 
alized objects, like galaxy groups and clusters, and the diffuse 
gas component which has low density and temperature. Therefore, 
the expected X-ray emission from the cosmic web can sensitively 
change when the parameters describing the ICM physics are varied. 
For these reasons, observational data on the soft X-ray background 
emitted by diffuse gas can be used as a probe of these physical pro- 
cesses , as suggested by different authors (se e, e.g.. lVoit & BrvarJ 
12001 d : lBrvan & VohfeOOltlXue & WJ2603I) . 

Thanks to the latest generation of X-ray satellites (Chandra 
and XMM-Newton), the observational picture about the cosmic X- 
ray background (hereafter XRB) in the soft (0.5-2 keV) energy 
band has become much clearer. Very deep pointed observations, 
like the 1-2 Ms Chand ra Deep Fields (CDFs) iGiacconi et all2002l 
Alexander et al. 2003) and the Lockman Hole XMM-Newton data 
IWorslev et al J 12004) . showed that the soft XRB is largely pro- 
duced by individual discrete sources, mostly AGN. Recent esti- 
mates of the percentage of the resolved c ontribution range be- 
tween 85 and 95 per cent (see, e.g jBauer et all20 04: Worsl ev et alJ 
120051) . These results give a conservative estimate of approximately 
(1.2 ±0.3) x 10~ 12 erg s cm~ 2 deg~ 2 for the intensity of the 
still unresolved soft X-ray background (see the discussion in Sec- 
tion ^. 1> : therefore this value can be used to set a stringent upper 
limit to the background from the diffuse IGM. 

The availability of these new observational constraints sug- 
gests to re-consider the problem of the soft X-ray emission from the 
cosmic web and its implications for th e ICM physics by updating 
an analysis made bv lCroft et alJ |2001>. The main goal of our paper 
is to check the consistency of our model with the observed limits 
by using the outputs of a cosmological hydrodynamic simulation 
( Borg ani et all2004T) , which, besides gravitational and gas dynam- 
ics, includes a treatment of the ICM with radiative cooli ng, star 
formation an d SN feedback . Our previous an alyses iBorgani et alJ 
l2004l lEttori et al Jl2004 iDiaferio et al JI2005I) showed that the re- 
sults of this simulation are in encouraging agreement with some 
of the most significant observed properties of clusters, such as the 
mass-temperature and the X-ray luminosity-temperature relations; 



however, the feedback mechanism does not appear to be efficient 
enough to avoid the production of overluminous groups and poor 
clusters. 

The paper is organized as follows. In Section|2|we describe the 
characteristics of our simulation and present the procedure used to 
create the set of X-ray maps. The general statistical properties of 
the XRB extracted from these maps are discussed in Sections|3]and 
|4|for the soft (0.5-2 keV) and hard (2-10 keV) X-ray bands, respec- 
tively. Section [5] is devoted to the comparison of our results with 
previous works. In Section|6|we review the observational estimates 
of the contribution of diffuse gas to the soft XRB and we compare 
them with the results of our mock maps in order to obtain qualita- 
tive constraints on the modeling of the physical processes included 
in our hydrodynamical simulation. We conclude in Section^ 



2 THE MAP CONSTRUCTION 

2.1 The hydrodynamical simulation 

The study of the properties of the diffuse X-ray emission from 
the large-scale structure of the universe and its comparison with 
observations requires the analysis of the total emission integrated 
over redshift. Therefore, we need to simulate the entire volume of 
the past light-cone seen by an observer located at z = 0. We con- 
struct the light-cone, which enables the production of simulated 
maps of the X-ray intensity, by using the outputs of a cosmological 
hydrodynamical simulation. Specifi cally, we use the results of the 
simulation by Borgani et al. (2004), which considers the concor- 
dance cosmological model, i.e a flat ACDM model dominated by 
the presence of the cosmological constant (£l m = 0.3, Q.\ = 0.7), 
with a Hubble constant Hq = h 100 km s _1 Mpc _1 and h = 0.7, 
and a baryon density £\ = 0.04. The initial conditions were set 
by the cold dark matter (CDM) power spectrum and were nor- 
malized by assuming Gg = 0.8. The run was carried out w ith the 
TreeSPH code GADGET-2 jSpringel et ai]|200ll: ISpringellEool 
and followed the evolution of 480 3 dark matter (DM) particles 
and as many gas particles from z = 49 to z = 0. The cubic box 
is 192/T 1 Mpc on a side, and the DM and gas particles have mass 
m DM = 4.62 x 1O 9 /i _1 M and m gas = 6.93 x 1O 8 /i _1 M , respec- 
tively. The Plummer-equivalent gravitational softening is £=7.5 h 
kpc at z = 0, fixed in physical units between z = 2 and z = 0, and 
fixed in comoving units at earlier times. 

The simulation, besides gravity and non-radiative hydro- 
dynamics, includes a treatment of the processes that influ- 
ence the physics of the ICM: star formation, by adopting a 
sub-resolution multiphase model for the interstellar medium 
ISpringel & Hernauistll2003l) . feedback from SNae with the effect 
of weak galactic outflows, radiative gas cooling and heating by a 
uniform, time-dependent, photoionizing UV background. This run 
produced one hundred outputs, equally spaced in the logarithm of 
the expansion factor, between z = 9 and z = 0. 

2.2 The map-making procedure 

In order to create mock maps of the X-ray emission of the structures 
wi thin the past light -cone, we follow the same technique adopted 
bv lCroft et alJ 1200 1) (see also similar applications for map s of the 
Sunyaev-Zerdovich eff ect, e.g. Ida Silva et all200lllSpringel et ail 
l200ltlWhite et al.l2002h . The method is based on the replication of 
the original box volume along the line of sight. Since we assume a 
flat cosmological model, we use the comoving coordinates for the 
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projection because the light trajectory is a straight line in this refer- 
ence frame. We build past light-cones which extend to z = 6. Even 
if this maximum redshift can appear too large to still yield a signif- 
icant signal in the X-ray band (this will be confirmed by our anal- 
ysis in Section[3}, this choice is done to obtain, at the same time, 
maps for the thermal and kinetic Sunyaev-Zel'dovich effects (not 
discussed in this paper), for which there is a non-negligible contri- 
bution from the high-redshift gas. The extension of the light-cone 
corresponds to a comoving distance of approximately 5,770 h 
Mpc, so we need to stack the simulation volume roughly 30 times. 
However, in order to obtain a better redshift sampling, rather than 
stacking individual boxes, we adopt the following procedure. We 
divide the simulated box at each output redshift into three equal 
slices along the line of sight (each of them with a depth of 64 
Mpc); for the stacking procedure, we choose the slice extracted 
from the simulation output that better matches the redshift of the 
central point of the slice. Our light-cones are thus built with 91 
slices extracted from 82 different snapshots. 

The necessity of avoiding the repetition of the same structures 
along the line of sight requires a randomization of the boxes used 
to build our maps: since our simulation assumes periodic bound- 
ary conditions, for each box entering the light-cone we combine 
a process of random recentering of the coordinates with a 50 per 
cent probability of reflecting each axis. The slices belonging to the 
same box undergo the same randomization process to avoid spa- 
tial discontinuities between them: this allows to retain the entire 
information on the structures within the box and strongly reduces 
the loss of power on larger scales. By varying the initial random 
seeds we also obtain different light-cone reconstructions: we use 
this technique to produce ten different realizations that we use to 
assess the statistical robustness of our results. 

In order to have maps covering a larger field of view, we repli- 
cate the boxes four times across the line of sight starting at comov- 
ing distances larger than half the light-cone extension (i.e. larger 
than about 2,900 h Mpc, corresponding to z > 1.4). The strategy 
to pile up boxes is shown in Fig. Q In this way we obtain maps 
3.78 deg on a side. Every map contains 8,192 pixel on a side: con- 
sequently our resolution is 1.66", roughly three times the resolution 
of the CDFs at the aimpoint. 

We now need to calculate the contribution, to the X-ray emis- 
sion, of every gas particle that lies within the light-cone volume. 
The X-ray luminosity of the ;-th particle in a given energy band 
[E1.E2], as measured at z = 0, is calculated as 

Lx,i = (/Jinpy 2 x e miPiA(Ti,Zi,E' l ,E' 2 ) , (1) 

where /j is the mean molecular weight in units of the proton mass 
m p , x e represents the ratio between the number density of free elec- 
trons and hydrogen nuclei (n e /nfj), rrij and p,- are the particle mass 
and density, respectively. Since only the gas with T > 10 5 K gives 
a significant contribution to the X-ray emission, we can safely as- 
sume full ionization of hydrogen and helium for all the particles: 
therefore the values of /j and x e depend only on the metallicity (no- 
tice that fj = 0.588 and x e = 1.158 for zero metallicity). The cool- 
in g function A is com puted by using the plasma emission model 
by Ravmond & Smith 1 1977) as a function of temperature 7) and 
metallicity Z; of the particle and depends on the energy interval 
[E[ , E^\ ; this interval is A"-corrected for redshift: E[ 2 — E\ 2 ( 1 + z) ■ 
As we will discuss in the following section, a sizeable contri- 
bution to the soft X-ray background comes from relatively warm 
particles, with temperature below 10 7 K. At these temperatures, 
the emissivity from metal lines becomes non negligible and, there- 
fore, one needs to account for their contribution. In its original im- 




Figure 1. Sketch of the configuration adopted to realize the light-cone. The 
observer is located at the position O at the centre of the lowest side of the 
first box. The past light-cone is obtained by stacking the comoving volumes 
of the simulation outputs at the corresponding redshift. In order to obtain 
a large field of view of size 3.78 2 deg 2 , starting at z = 1.4 we use four 
replications of the box at the same redshift. The red lines show the volume 
in the light cone corresponding to the field of view. 



plementation, GADGET-2 includes a prescription to generate met- 
als from SN explosions. This model assumes that only SN-II con- 
tribute to the chemical enrichment, under the assumption of instan- 
taneous recycling, i.e. metals are released instantaneously when 
new stars form and the effect of stellar life-times is neglected (see, 
e.g., |jornatore et al. 2004, for a more detailed implementation of 
chemical enrichment in GADGET-2). In Eq. Q we adopt the value 
of metallicity Z; yielded by this prescription to compute the lumi- 
nosity of each gas particle. 

We then calculate the contribution of the i-th gas particle to 
the X-ray intensity Ixj as 

Ix,i=L X j/(4%d L (z) 2 A) , (2) 

where di(z) is the luminosity distance and A is the angular area 
covered by a pixel. This quantity is then distributed over the pixels 
by using an SPH smoothing kernel given by 
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Figure 2. Maps of the soft (0.5-2 keV) X-ray intensity obtained by considering all the gas particles (IGM, left panel), and only the gas particles with 
temperature in the range 10 5 < T < 10 7 K (WHIM, right panel). The maps are 3.78 deg on a side and the pixel size is (1.66") , These two maps refer to the 
same realization of the past light-cone. 



( l-6x 2 +6x 3 , (Kx<0.5 
W(x)oc) 2(1 -x) 3 , 0.5<x<l (3) 
{ 0, x> 1. 

In the previous expressions x = A9/oc;, where AO is the angular 
distance between the pixel centre and the projected particle posi- 
tion and a,- is the angle subtended by the particle smoothing length 
provided by the hydrodynamical code. In order to conserve the total 
intensity emitted by each particle we normalize to unity the sum of 
the weights W over all "touched" pixels. Finally, the intensity for a 
given pixel is obtained by summing over all the particles inside the 
light-cone. 

To avoid spurious effects in the computation of the X- 
ray intensity, we exclude the particles having a mean elec- 
tron density n e > . 26/i 2 c m~ 3 . According to the model of 
SDi'insel & Hernauist 1 2003J, these particles are assumed to be 
composed by a hot ionized phase and a cold neutral phase, whose 
relative amounts depend on the local temperature and density. Since 
these particles are meant to account for the multi-phase nature of 
the interstellar medium, we correctly exclude them from the com- 
putation of the X-ray emissivity. 



3 THE PROPERTIES OF THE SOFT X-RAY 
BACKGROUND 

Examples of intensity maps obtained by adopting the method de- 
scribed above are displayed in Fig. |2|f or the emission in the soft 
(0.5-2 keV) band. Following lCroft et alJi200lT) . we distinguish two 
X-ray contributions: the contribution from the intergalactic medium 
(IGM, left panel), i.e. from all the gas particles, and the contribution 
from the warm-hot intergalactic medium (WHIM, right panel), i.e. 
from the gas particles having a temperature between 10 5 and 10 7 K. 
The IGM map is dominated by several extended bright galaxy clus- 
ters (with flux in the range 10~ 12 — 10~ n erg s _1 cm -2 ) that give 
significant contributions to the total flux; there is also a large num- 
ber of smaller structures, corresponding to nearby galaxy groups or 



distant fainter clusters. These smaller objects give the dominating 
contribution in the WHIM maps, because the WHIM does not con- 
tain the hottest particles; on the contrary, in these maps, the bright- 
est galaxy clusters are much less prominent. Note that, because of 
projection effects, neither map shows evident signatures of emis- 
sion from filamentary structures. 

The distribution of the values assumed by the intensity (con- 
sidering the 1.66"pixels) is shown in Fig.|5|for both the IGM and 
the WHIM. The thin lines refer to each of the ten light-cone real- 
izations, while the thick solid line shows their average. A small dis- 
persion between the different realization is evident: this is mainly 
due to the inclusion, in the maps, of a larger or smaller number of 
galaxy clusters at relatively low redshifts. In general the two av- 
eraged distributions for / = loglx (erg s _1 cm~ 2 deg -2 ) are very 
close to a Gaussian distribution. In particular we find / = —12.81 
with a corresponding r.m.s. of 1.05 for the IGM, and / = —13.01 
with a corresponding r.m.s. of 0.94 for the WHIM; the skewness is 
0.30 and 0.25, for the IGM and the WHIM respectively. Of course 
the distributions differ for larger fluxes: this is due to the fact that 
the flux of the brightest structures comes from gas at high temper- 
ature. 

The values of the intensity averaged over the ten map realiza- 
tions is reported in TableQfor both the IGM and the WHIM. The 
quoted errors are the dispersions computed in fields of 1 deg 2 to 
al low a direct comparison with the corresponding results reported 
in lCroft et ail 1200 ll) . We find that the mean intensity for the IGM 
is about 4 x 10~ 12 e rg s cm~ 2 deg -2 , a factor 1.8 larger than 
in lCroft et alj 1200 ll) . The corresponding dispersion is also larger: 
about 50 per cent compared to about 20 per cent in ICroft et aU 
1200 ll) . Again, we checked that the high map-to-map spread is due 
to the presence of some bright galaxy clusters very close to the ob- 
server and giving a strong contribution to the total emission: this 
effect is also evident in the pixel distributions of individual maps 
shown in Fig.|3| By looking at the WHIM contribution, we find that 
its mean intensity is about 40 per cent of the total bremsstrahlung 
emission in the soft band (see Table0. Again this value is larger 
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Figure 3. Distribution of pixel values of the intensity in the soft X-ray band. 
Upper and lower panels refer to the IGM and the WHIM, respectively. The 
thin lines show the results often different realizations; the thick solid line 
is the corresponding average. 



Table 1. Average value of the intensity Ix in the two different X-ray bands 
for the IGM and the WHIM. The average is computed over ten different 
map realizations; the quoted errors are the r.m.s. in fields of 1 deg 2 . 



Energy band 


IGM (total) 


WHIM 






(erg s _1 cm -2 deg -2 ) 


(erg s cm -2 


deg 2 ) 


Soft (0.5 - 2 keV) 


(4.06 ±2.19) x 10~ 12 


(1.68 ±0.62) > 


;10- 12 


Hard (2 - 10 keV) 


(1.01 ±1.53) x 10~ 12 


(2.92 ±2.46) > 


; 10- 14 



than in the analysis of ICroft et alJ 1200 ll) . who reported a value 
of 4.15 x 10~ 13 erg s _1 cm -2 deg -2 , corresponding to about 20 
per cent of the total intensity of the IGM. However, we emphasize 
that we cannot make a direct c ompari son between these results and 
those presented in lCroft et alJ 1200 ll) because the simulation ana- 
lyzed in that paper considers a much smaller box (50h~ l Mpc on a 
side), assumes a slightly different cosmological model (a flat uni- 
verse with Q. m = 0.4 and a primordial spectral index n = 0.95) and, 
most importantly, the included physical processes are not the same 
(see the discussion in[5j. 

To determine which sources contribute most to the soft X-ray 
emission, in Fig.[4]we show how the mean intensity varies when we 
include pixels with fluxes above a given threshold. The results show 
that the bremsstrahlung emission is dominated by very bright ob- 
jects: pixels having a surface brightness larger than 10~ 10 erg s" 1 
cm -2 deg -2 (mainly corresponding to galaxy clusters and groups 
as shown in Section l6"2l contribute to more than 50 per cent of the 
total value, for both the IGM and the WHIM. 

In order to have an indi cation of the typic al distances of the 
soft X-ray sources, we follow ICroft et alJ 1200 ll) and compute, for 
each pixel, the mean redshift of the particles contributing to the 
XRB weighted by their fluxes. The resulting map for the same 
IGM realization shown in the left panel of Fig.|2|is displayed in 
Fig.[5](left panel). Almost all the bright clusters are at low redshift 
(z < 0.15), while the emission from the faintest structures mainly 
comes from redshifts around unity. This is confirmed by the dis- 
tributions of the flux-weighted redshifts (central panel of Fig. [5}> 



Figure 4. The value of the mean intensity (averaged over ten different map 
realizations) as a function of the lower limit of the pixel surface brightness. 
The solid and dashed lines refer to the IGM and the WHIM, respectively. 

computed from the complete set of ten realizations: the shapes of 
the distributions for the IGM and the WHIM are very similar, with 
a peak at z ~ 0.7. Notice that few realizations have an excess of 
contribution from sources at relatively low redshift (z < 0.15): this 
is the origin of the large spread of the mean intensity previously 
discussed. The analysis of the maps also clearly shows how the pro- 
jection of nearby and distant objects makes it very difficult to see 
imprints of the filamentary network present in the large-scale struc- 
ture of the universe. Finally, in the right panel of Fig.|5|we show the 
integrated soft X-ray intensity as a function of the limiting redshift. 
The curves for both the IGM and the WHIM reach half the total 
value at z — 0.3, the 90 per cent level is reached around z — 0.9 for 
the IGM and at z ~ 0.8 for the WHIM. We can also notice that the 
sources at z > 2 contribute less than 1 per cent of the soft X-ray 
flux. 



3.1 The autocorrelation function 

We analyze the clustering properties of the IGM in the soft X-ray 
band by estimating the autocorrelation function in our simulated 
maps. Following Croft et al. 1 2001), we calculate the contrast 8 (it) 
defined as 

5(x) = [/ x (x)/7]-l , (4) 

where Ix (x) is the surface brightness at the position x and / is the 
average of each map. The contrast is then used to compute the an- 
gular correlation function as 

w(6) = (8(x)8(x+9)> . (5) 

The results are shown in Fig. [6] where we plot the average 
of the angular correlation functions calculated by considering the 
ten different realizations of the light cone; error bars are the error 
on the mean obtained from the scatter of the maps. We notice that 
the mean function is always positive up to 6 ~ 30', although some 
of the realizations show negative values for > 8'; the smallest 
scale at which w(0) is compatible with a null value is 20'. From the 
plot it is also evident a slow change of the logarithmic slope. This 
is confirmed by our analysis: by fitting our mean function under 
the assumption of a power-law relation, w(9) °< 9 a , we find a ~ 
-0.9 and a w -2.6 in the intervals 0.1' < < 2' and 2' < < 30', 
respectively. 
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Figure 5. Left panel: the map of the flux-weighted redshift, computed for the IGM in the soft X-ray band. The plot refers to the same light-cone realization 
shown in Fig. [3] Central panel: as in Fig.|5] but for the distribution of the flux-weighted redshift in the soft X-ray band. Right panel: the integrated soft X-ray 
intensity (normalized to the total value) as a function of redshift. The solid and dashed lines refer to the IGM and the WHIM, respectively. The values plotted 
are the mean of ten map realizations. 
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Figure 6. The angular correlation function of the IGM in the soft (0.5-2 
keV) X-ray band. The solid line is the average over ten different maps 3.78 
deg on a side; error bars are the error on the mean calculated from the scatter 
of the maps. 

Our results can be directly compared to those obtained by 
ICroftetalJiioOlll . In general the autocorrelation function extracted 
from our maps is larger by a factor of about 6. Moreover we find 
that at large scales the decrease of their w(0) is slower (they found 
a = —1.4). These discrepancies can be easily explained by con- 
sidering the differences between the two simulations . First, and 
most i mportantly, the size of the simulation box used bv lCroft et alJ 
1200 ll) is about four times smaller than our box size: the resulting 
lack of large-scale power can lead to a significant underestimate of 
the clustering strength. Second, the physical processes included in 
our simulation produce higher values of the soft X-ray signal: as a 
consequence, the contrast S also tends to be amplified, thus origi- 
nating a larger autocorrelation function. 



4 THE CONTRIBUTION IN THE HARD X-RAY BAND 

It is well known that the hard (2-10 keV) X-ray band is strongly 
dominated by the emission from AGN. However, a small portion 
of gas, mainly located within the richest clusters, is at high tem- 



perature and also gives some contribution in this band. The maps 
presented in Fig.Qshow the hard band intensity for the IGM and the 
WHIM (left and right panels, respectively) for the same light-cone 
realization shown in Fig.|2| We note that, even if in some positions, 
corresponding to the hottest clusters, the intensity can reach val- 
ues similar to those obtained in the soft X-ray maps (10~ 9 erg s" 1 
cm" 2 deg -2 ), outside virialized objects the hard X-ray maps are 
much fainter. In the map for the WHIM the signal is strongly re- 
duced and limited to the cluster atmospheres, with no evident emis- 
sion from the diffuse gas: this makes the cosmic web not observ- 
able in this band because of its very low signal-to-noise ratio and 
projection effects. However, the high contrast between clusters and 
diffuse regions in the hard band maps suggests an empirical way to 
identify the positions of galaxy clusters in our simulated maps. This 
method will be used in Section|6|to separate the emission coming 
from extended objects from that produced by diffuse gas. 

The small contribution of the IGM and the WHIM to the hard 
XRB is confirmed by the average total intensity reported in TableQ 
for both components: we find about 1 x 10~ 12 erg s" 1 cm" 2 deg" 2 
and 3 x 1(T 14 erg s" 1 cm" 2 deg" 2 , for the IGM and the WHIM, 
respectively. These values account for approximately 25 per cent 
and 2 per cent of the corresponding mean intensity in the soft band. 
Note that, while for th e WHIM our flux is a factor of two larger 
than ICroft et alJ 120011) result, for the IGM there is a substantial 
ag reement. In fact, th e lack of large-scale power in the simulation 
bv ICroft etaP 1200 ll) . due to the small box size, is compensated 
by the higher dynamical evolution associated to the higher value 
of Q. m , for a fixed power spectrum normalization Os assumed in 
both analyses. Moreover, the effect of feedback on the hard X-ray 
emission by rich clusters is expected to be less relevant. We have 
to consider, however, that the field-to-field variance is quite large 
for the hard X-ray band, because most of the signal comes from 
nearby bright clusters: in fact, the hard X-ray spectrum decreases 
exponentially and, even for hot objects, the bulk of the signal is 
rapidly shifted into the soft band at increasing redshifts. 

A way to discriminate the AGN signal from that coming 
from clusters and the cosmic web is based on the energy flux 
hardness ratio (HR), here defined as the ratio between the inten- 
sities in the hard and soft bands at each pixel position: HR = 
I x [2 - 10keV]// x [0.5 - 2keV]. The emission from AGN is in fact 
expected to have HR > 0.5, while much lower values are common 
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Figure 7. Maps of the hard (2-10 keV) X-ray intensity for the IGM and the WHIM (left and right panels, respectively). The map is 3.78 deg on a side and the 
pixel is 1.66" on a side. These maps have been obtained from the same light-cone realization shown in Fig.l2l 



for galaxy clusters and diffuse gas. The map for the HR obtained 
from the same IGM realization displayed in Figs.[2|and0is shown 
in the left panel of Fig. [8] in the plot it is easy to recognize the posi- 
tion of the richest galaxy clusters, while the connecting filamentary 
structure is almost completely absent. The distributions of the HR 
values for the ten realizations and their average are shown in the 
right panels, for both the IGM (up) and the WHIM (down). In the 
former case the median value is around HR = 0.03, and less than 
a few per cent of the pixels have HR > 0.5; in the latter case the 
median is smaller than 0.01 and no pixel has HR > 0.1. 



5 COMPARISON WITH PREVIOUS WORKS 

In this paper we study the properties of the X-ray background 
from the large-scale structure of the universe. In particular we fol- 
lowed the same method to build light-cone realizations applied by 
ICroft et al. (2001), performing then similar statistical tests. From 
a comparison between our and their analysis, we found that in our 
soft X-ray maps the mean intensities of the IGM and the WHIM are 
larger than the intensities obtained bv lCroft et alj 1200 ll) by a fac- 
tor of 1.8 and 4, respectively. The autocorrelation function of the 
signal is also larger by a factor of about 6. As for the hard X- ray 
band, we find a substantial agreement with lCroft et alj 1200 ll) for 
the IGM contribution, while our contribution from the WHIM is a 
factor of 2 higher. 

These differences can be understood in terms of the different 
characteristics of the two hydrodynamic simulations on which the 
analyses are based. The main reason for our higher values is that we 
followed the evolution of a much larger (almost 60 times) volume. 
This allows us to obtain a better representation of very massive 
clusters that contribute to most of the emission (see the discussion 
in Section [3j. The larger sampled volume is also the main moti- 
vation of our higher autocorrelation function, which includes the 
contribution of longer wavelength modes. 

However, an important difference between the two hydrody- 
namic simulations is also related to the treatment of the physical 



processes affecting the e volutio n of the baryonic component. The 
simulation of Croft et al. 1 2001 ) did not consider the possible mul- 
tiphase nature of the starforming gas. For this reason they had to 
resort a posteriori to a correction of their results to account for it. 
On the contrary, the GADGET-2 code used for our simulation fol- 
lows the different physical processes (radiative cooling, star forma- 
tion and supernova feedback) in a self-consistent way. In particular 
the inclusion of a phenomenological description of galactic winds 
provides a much more efficient energy feedback than that imple- 
mented in the simulation by Croft et al. 1 2001 ). Moreover the code 
directly considers the effect of a photoionizing time-dependent UV 
background. Globally this gas treatment produces a larger fraction 
of WHIM. This explains why larger differences are found for the 
contributions of the WHIM, while less significant differences are 
found between the signals in the hard X-ray band, which are mostly 
originated by the hottest objects. We also assumed a cosmolog- 
ical model which is closer to that sugg ested by the mos t recent 
observational data with £2 m = 0.3, while C roft et all200ll adopted 
£l m = 0.4: thus we have a higher contribution from baryons at high 
redshifts due to the slower dynamical evolution of the structures. 
The apparent agreement o f our re sult for the hard band emission of 
the IGM with lCroft et alj 1200 ll) is due to a combination of all of 
these factors (see the discussion in Section ^. 

Finally, a simila r result to the one o f lCroft et all 1200 ll) has 
been obtained also by Bryan & Voit ( 2001) for the soft X-ray emis- 
sion of the IGM from their AMR simulation which includes a 
simple feedback model. Again the difference with respect to the 
present work is maninly due to the relatively small volume sam- 
pled by their simulation. 



6 THE CONTRIBUTION OF DIFFUSE GAS TO THE 
SOFT X-RAY BACKGROUND 

The goal of this section is to extract, from the soft X-ray maps dis- 
cussed above, the contribution produced by the cosmic web, i.e. the 
baryons which are not included in the groups and clusters of galax- 
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Figure 8. Left panel: the map of the hardness ratio (hard band to soft band) HR. This map is obtained from the same realization of the IGM light-cone shown 
in the left panels of Figs.l2landl71 Right panels: distribution of HR for the IGM (up) and for the WHIM (down). The thin lines show the results for ten different 
realizations; the thick solid line is the corresponding average. 



ies currently identified in deep X-ray surveys. In order to check 
whether the modelling of the physical processes treated by our hy- 
drodynamical simulation (which also influence the low-density gas 
regions) is reliable, we compare our results with observational esti- 
mates (mainly upper limits). We have to notice that in the literature, 
as reviewed in the following subsection, there is no robust determi- 
nation of the contribution of diffuse gas to the soft XRB: the value 
that we will use later for comparison is based, to some extent, on 
interpretation of the existing data. 

6.1 Observational estimates 

In recent years, thanks to the availability of new X-ray satel- 
lites like Chandra and XMM/Newton, there have been many ef- 
forts to measure the tot al XRB and to resolve it in its differ- 
ent components (see, e.g. iGiacconi et all2002LlMoretti et alJ2002L 
lAlexander et all200llworslev et al 120041: iBauer et all2004 and 
references therein). In fact the possibility of recognizing the dis- 
crete sources which give the main contribution opens new windows 
to study the connection between AGN and galaxy formation and, 
most importantly for this work, constrains the possible signals from 
the diffuse gas. 

At present, the most reliable estimate of the total XRB in the 
soft (0.5-2 keV) band has been obtained bv lWorslev et alJ 12005T) . 
By assuming the shape of the XRB in the [1-8 keV] band deter- 
mined bv iDe Luca & Molendil 12004) from BeppoSAX data, and 
by taking into accou nt the steepening of t he spectrum below 1 
keV as found bv lRoberts & Warwick! ^OOlh . lWorslev et alj<2005l) 
estimate an XRB[0.5 - 2keV] = (8.12 ±0.23) x 10~ 12 erg s" 1 
cm~ 2 deg~ 2 . Note that this value is larger than the one adopted by 
IBauer et alJ J2004I) . who found XRB [0.5 - 2keVl = (7.52 ±0.35) x 
10~ 12 er g s -1 cm~ 2 deg~ 2 . The discrepancy is originated by the 
fact that IBauer et alJ 120041) extrapolate the hard spectrum in the 
soft band by assuming a slope of —1.4: of course, this does not 
take into account the steepening due to the non-AGN contribution, 
which, by the way, is the one we are most interested in. This to- 



tal value includes any contribution, namely the Local Hot Bub- 
ble, the Galactic halo, unresolved galactic stars, and most impor- 
tantly, the extragalactic components (both bright and faint) from 
point sources (AGN), clusters, groups and the truly diffuse gas. 
The resolved (to date) XRB in the soft band prod uced by these 
extrag alactic sources has been recently estimated bv lWorslev et alJ 
(2005), who found (6.9 ±0.2) x 10~ 12 erg s" 1 cm~ 2 deg~ 2 . Simi- 
lar res ults for the resolved component are presented bv lBauer et alJ 
120041) : (6.7 ±0.3) x 10~ 12 ergs' 1 cm' 2 deg~ 2 . Therefore, by as- 
suming the values found bv lWorslev etafl <2005h for the total XRB 
in the soft band, the emission which is still unresolved is approx- 
imately (1.2 ±0.3) x 10~ 12 erg s _1 cm~ 2 deg~ 2 . This value rep- 
resents an upper limit to the possible contribution from the diffuse 
gas, because it still includes the unknown Galactic and Local Hot 
Bubble contributions. Notice that the latter components probably 
dominate the XRB below 0.5 keV, making highly uncertain the es- 
timate of the extragalactic contribution in the lowest X-ray energy 
band, where the signal from the cosmic web is expected to be even 
more significant. 



6.2 Separating the diffuse soft X-ray emission from galaxy 
groups and clusters in the simulated maps 

By comparing naively the estimates of the IGM soft X-ray con- 
tribution reported in TableQwith the observational data discussed 
above, we would notice that our simulated XRB accounts for about 
50 per cent of the observed value of the XRB, but it is much larger 
(by a factor 4) than the observed upper limit obtained by removing 
the discrete sources. Even considering the signal from the WHIM 
only, i.e. from gas having temperature between 10 5 and 10 7 K, the 
signal in our maps would be too large. However, this comparison 
is quite misleading, because in the results from the simulation we 
still include the portion of the warm gas located in groups or in 
the external atmospheres of galaxy clusters; on the contrary, this 
gas is not present in the observed value 1.2 ±0.3 x 10~ 12 erg s" 1 
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cm deg . For this reason we need a more detailed analysis of 
our maps, where these objects need to be excluded. 

We already know, from the previous analy sis of our original 
hydrodynamical simulation (see, e.g., Borgani et al. 2004), that the 
modeling of the included physical processes produces overlumi- 
nous groups and poor clusters of galaxies. As discussed in that 
paper, additional sources of energy (like feedback from SN-Ia or 
AGN) would help to solve this discrepancy. This excess is located 
in high-density regions, whereas we are interested in the soft X-ray 
emission from the cosmic web: therefore, we need to identify and 
remove the X-ray extended sources, which correspond to galaxy 
clusters and groups, from the total intensity of our maps. 

The proper procedure would be to identify and remove all the 
clusters and groups from mock observations of our simulated vol- 
umes. In practice, this is unfeasible, because it would require sim- 
ulating a full wide-area survey with all the related observational 
aspects, which is far beyond the aim of this paper. Therefore, we 
will remove clusters and groups on the basis of their total flux as it 
appears on our flux maps. To establish a definition of a cluster or a 
group to be applied to our bidimensional maps, we use an identifi- 
cation criterion based on the surface brightness value. We identify 
as an X-ray extended source (group or cluster) all the connected 
regions of pixel above a surface brightness threshold. 

To accomplish this, we use the maps in the hard (2-10 keV) 
band because the emission from the diffuse gas in this band is much 
lower than in the soft band and therefore it is easier to identify the 
regions corresponding to virialized haloes (see Fig.Q. Then, we 
choose the value of the surface brightness threshold by inspecting 
the faintest groups identified in the deepest X-ray observations to 
date. Specificall y, we consider a gro up (CDFS-594) identified in 
the CDF South iGiacconi et all2002l) at z ~ 0.7 - 0.8. Its emission 
is detected with a high signal-to-noise ratio over a circular region of 
approximately 2700 a rcsec 2 . After fi tting its spectrum with a mekal 
model (see, e.g-. lLiedahl et alJl995l and references therein) we ob- 
tain a temperature of about 2. 1 keV. Therefore, the corresponding 
average surface brightness in the hard band is 3.06 x 10~ 12 erg s — 1 
cm~ 2 deg -2 , which will be used as the reference value in the fol- 
lowing analysis 1 . Notice that this average value corresponds to the 
central part of the group (typically about 3 core radii), which is the 
region where the emission can be realistically detected. However, 
some emission associated to the X-ray halo is expected to provide 
sub-threshold contribution up to the virial radius. Since we want to 
remove all the X-ray emission associated to a given halo, in the fol- 
lowing analysis we will also consider smaller, up to 8 times, surface 
brightness thresholds. We choose these thresholds, because the sur- 
face brightness of an isothermal gas in the central region (about 3 
core radii) is roughly 8 times larger than at the virial radius. Finally, 
we consider a threshold for the minimum angular size of the region 
that can be identified as a group or a cluster. For this purpose we 
choose a reference physical length of 100 kpc, that corresponds to 
a minimum angular size of ~10 arcsec (at z ~ 1) for the cosmologi- 
cal model assumed in our simulation: this gives a minimum angular 
surface of ~3 10 arcsec 2 . This quantity roughly corresponds to the 
minimum size of the central detectable regions of X-ray groups (see 
IWillis et alj2005l) . This last criterion has a small impact on the final 
results (i.e. the expected average surface brightness of the cosmic 
web changes of few percent when doubling the minimum angular 
size of the haloes). 



Similar results are obtained by considering a different galaxy group 
(CDFS-645), again identified in the CDF South exposure. 



6.3 The soft X-ray emission from the cosmic web in the 
simulated maps 

By applying the previously discussed thresholds of the surface 
brightness and size on the maps produced in the hard (2-10 keV) X- 
ray band, we create a catalogue of connected regions corresponding 
to the extended objects like groups and clusters. We use them as a 
mask on the soft (0.5-2 keV) X-ray maps to remove their emission 
and finally obtain the diffuse contribution. By averaging over our 
ten different realizations, we find that the signal from the diffuse 
gas is 1.57 x 10~ 12 erg s _1 cm~ 2 deg~ 2 which is comparable to 
the observed upper limit. 69 per cent of the diffuse emission comes 
from gas with temperature in the range 10 5 — 10 7 K, whereas the to- 
tal contribution of clusters and groups to the XRB is 2.49 x 10~ 12 
erg s _1 cm -2 deg -2 ; this value is larger than the total contribu- 
tion from cluste rs and grou ps estimated from the observed number 
counts iRosati et alJl2002l see also the corresponding discussion 
at the end of this section). As noted above, the adopted surface 
brightness threshold can be considered as an average of the group 
emission at intermediate redshifts, and then it could be too high to 
account for all clusters' and groups' emission. In fact we find that 
31 per cent of the 1.57 x 10~ 12 erg s _1 cm~ 2 deg~ 2 emission is 
still produced by gas having a temperature larger than 10 7 K. This 
is further confirmed by a visual inspection of the masked maps: 
the external regions of nearby clusters are not excluded. This is an 
indication that in general the surface brightness threshold defining 
high density regions associated to groups and clusters of galaxies 
is probably lower than the adopted value. For this reason, we re- 
peat the previous analysis by using thresholds which are smaller 
by factor of 2, 4 and 8. The results, reported in Table [2] show a 
reduction of the contribution of the diffuse gas by a factor of 2 
when the surface brightness cut is 8 times smaller; the correspond- 
ing fraction of the signal from the WHIM increases up to 82 per 
cent. In general we can conclude that the soft X-ray emission from 
the diffuse gas computed from our simulation is consistent with the 
upper limit coming from observational data. The feedback model 
we are assuming is efficient in the WHIM and/or in regions with 
intermediate density, although it can show problems in the deepest 
potential wells. 

Finally, as a self-consistency check, we compute the logA'- 
logS function (always in the soft X-ray band) for our connected 
regions. The results are shown in Fig.|9|for the four different thresh- 
olds reported in Table [2] Again, we notice the excess with respect 
to the observati onal data, represented by a combination of data 
made bv lRosati et all 1 2 002) and coming from the Rosat Deep Clus- 
ter Survey iRosati et alJll998h . EM SS iRosati et alj[l 995). BCS 
lEbeling et alJll998t) and REFLEX jBoehringer et alj|200ll) . This 
excess is produced by the overluminous objects which are present 
in our simulation at all relevant redshifts: for this reason, number 
counts are overestimated over the entire flux range. Finally, the ap- 
parent flattening at low fluxes (S ~ 10~ 15 erg s _1 cm~ 2 ) is orig- 
inated by both resolution effects and the surface brightness limit 
assumed for the identification of the extended sources. 



7 CONCLUSIONS 

In this paper we have used a cosmological hydrodynamical sim- 
ulation of a concordance ACDM m odel to discu ss the properties 
of the diffuse XRB. The simulation l Borgani et al. 2004) includes 
many relevant physical processes affecting the gas component. In- 
deed the numerical treatment accounts for a time-dependent, pho- 
toionizing UV uniform background, radiative cooling processes 
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Table 2. Contribution to the soft (0.5-2 keV) XRB of galaxy clusters and groups (Column 2) and of the diffuse gas (Column 3) as a function of the adopted 
threshold of the surface brightness (Column 1). We also report the percentage of the diffuse contribution coming from the WHIM (Column 4). The quoted 
eiTors are the r.m.s. in fields of 1 deg 2 . All values are in units of 10"' 2 erg s"' cm -2 deg" 2 . 



Soft (0.5-2 keV) X-ray background 
(10" 12 ergs" 1 cm" 2 deg" 2 ) 
SURFACE BRIGHTNESS THRESHOLD CONTRIBUTIONS FROM 



IN THE HARD (2-10 KEV) X-RAY BAND 


CLUSTERS AND GROUPS 


Diffuse gas 


%WHIM 


3.06 


2.49 it 2.04 


1.57 ±0.27 


(69 ±13)% 


1.53 


2.78±2.10 


1.28±0.19 


(73 ±12)% 


0.77 


3.05±2.14 


1.01 ±0.13 


(78 ±11)% 


0.38 


3.28±2.17 


0.78 ±0.09 


(82 ±11)% 


OBSERVATIONAL UPPER LIMIT 




1.2 ±0.3 






IQ-15 10 -14 10" 13 10-' 2 10-' 1 

S[o,5-2] (erg s-' cm- 2 ) 



Figure 9. Number of objects as a function of the limiting flux in the 
soft (0.5-2 keV) X-ray band. The black solid lines are the results obtained 
by considering the connected regions with four different surface bright- 
ness cuts (see Table[2]for the adopted values). For comparison, we show 

^ I i I j ■ 

(red dashed line) the observational data obtained by Rosati et al. 1 2002 ) by 
comb ining the results from the Rosat Deep Cluster Surve y i Rosati et al. 
Il998t). EMSS iRosati et alll995l) . BCS lEbeling et alll998h and REFLEX 
i Boehringer et all200ll) . 

within an optically thin gas of hydrogen and helium in collisional 
ionization equilibrium, star formation events, and feedback pro- 
cesses from both SN-II and galactic winds. Our previous analy- 
ses iBorgani et al 120041: lEttori et all2004l) showed that the results 
of this simulation are consistent with several observational X-ray 
properties of galaxy clusters. However, there is also a number of 
discrepancies that remain unaccounted for: in particular the cluster 
X-ray luminosity-temperature relation extracted from our simula- 
tion appears too high, i.e. there are overluminous groups and small 
clusters. 

As suggested by different authors (see , e.g., IVoit & Brvanl 
1200 id: iBrvan & Voitll200lt IXue & WulEooH) . an alternative way 
to constrain the model describing the thermal properties of baryons 
and their cosmic history is related to the soft X-ray emission from 
diffuse gas. Indeed, a model of heating and cooling of the inter- 
galactic medium makes predictions in terms of diffuse emission 
from filaments and unresolved structures. These predictions can be 
compared with the observational measurements of the XRB, which 
are now available (in form of upper limits) thanks to the deep Chan- 
dra data. 



In order to test t he model i n cluded in our simulation, we 
follow the method of ICroft et alJ 1200 ll) and, starting from the 
simulation snapshots, we constructed a set of ten different two- 
dimensional maps, of size (3.78deg) 2 of the past light -cone back 
to z = 6. The present analysis extends Icroft et all 1200 ll) previous 
work, thanks to a much larger volume sampling, to a more realistic 
representation of the physical processes affecting the baryon his- 
tory and to the comparison with the observed XRB data after the 
Chandra era. 

The main results obtained by using our set of maps are as fol- 
lows: 

• The mean intensity of the IGM in the soft (0.5-2 keV) X-ray 
band is about 4.1 x 10" 12 erg s" 1 cm" 2 deg" 2 ; when considering 
the WHIM (defined as gas with temperature between 10 5 and 10 7 
K) the mean intensity reduces to about 1.7 x 10" 12 erg s" 1 cm" 2 
deg" 2 . The distribution of the intensity in the maps is very similar 
to a lognormal. 90 per cent of the signal comes from structures at 
z<0.9. 

• As expected, in the hard (2-10 keV) X-ray band the total mean 
intensity is smaller (by a factor 4) than in the soft one and becomes 
almost negligible when considering the WHIM (about 3 x 10" 14 
erg s _1 cm -2 deg -2 ). The hardness ratio in the maps has a dis- 
tribution which peaks at low values, with a median close to 0.03, 
which enables an easy discrimination of the AGN signal. 

• We obtain an estimate of (0.8 — 1.6) x 10" 12 erg s _1 cm -2 
deg" 2 for the soft X-ray emission from the diffuse gas after re- 
moving the regions corresponding to extended objects, like galaxy 
groups and clusters, from the maps. Our result is consistent with 
the present upper limit coming from observational data: (1.2 ± 
0.3) x 10" 12 erg s" 1 cm" 2 deg" 2 . This value, recently obtained by 
Worsl ev" st alJ i2005). has been estimated by removing the contri- 
bution of resolved sources (AGN, groups and clusters of galaxies) 
from the total XRB. 

As a conclusion, our results show that the physical processes 
discussed here are consistent with existing constraints on the X- 
ray properties of the warm-hot baryons. However, if in the near 
future the measurements of the possible contribution from Galactic 
and local structures (still included in the observational upper limits) 
will reduce the maximum allowed emission from the diffuse gas by 
a factor of 2 or more, our predictions will start to conflict with the 
data. In this case it will be required a more efficient mechanism 
of feedback acting on the gas with intermediate temperatures and 
densities: in fact the presence of a stronger feedback would place 



11 



the gas on a higher adiabat, thereby preventing it from reaching 
high densities and, on turn, suppressing the diffuse emission. 

Our analysis confirms that the comparison between observa- 
tional constraints on the diffuse emission in the soft X-ray band and 
results from cosmological hydrodynamical simulations enables us 
to gain information on the thermodynamical history of the diffuse 
baryons. In the future it will be interesting to study some detailed 
simulated observations of the diffuse gas to check its detectability 
not only through bremsstrahlung emission, but also via emission 
and/or absorption of oxygen lines (Ovi, Ovil and Ovill): real ob- 
servations will probably be provided by some dedicated wide-field 
high-resolution spectroscopic experiments now in project. It will 
also be relevant to compare simulated results with observations per- 
formed by the next generation of X-ray spectroscopic satellites, like 
Constellation-X/XEUS or NeXT, provided that their field of view is 
large enough to cover the filamentary structure of the cosmic web. 
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